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Abstract 

We consider a description of the deuteron based on meson exchange potentials. A key feature is the 
inclusion of the I = S = two-pion intermediate state ('cr(600)') as a significant component of the inter- 
nucleon potential energy. In this approach, deuteron binding is seen to be predominantly a consequence of 
er(600) and w(783) exchange, with a secondary role played by p(77Q). We explore sensitivity of two-nucleon 
binding to changes in the potential and thereby obtain an anthropic constraint — that the deuteron unbinds 
for a modest decrease (about 6%) in the attractive cr(600) potential. 
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I. INTRODUCTION 



Imagine a sequence of worlds in which the light quark masses are continuously varied away from 
their physical values. It has been argued that heavy nuclei will disassociate for a 64% increase in 
the sum m u + vn^ [1]. One can anticipate that before this happens, the deuteron will become 
unbound. This is because the deuteron binding energy is less than the average binding energy per 
nucleon. Here, we study this problem quantitatively and provide details about how the deuteron 
would respond to variations in its potential induced by changes in m u + [2]. 

The deuteron is an example of a hadronic molecule. Such a two-hadron composite, bound by the 
strong interactions and whose mass lies slightly below the associated two-hadron threshold, should 
be describable with a nonrelativistic Schrodinger equation. We can picture the two hadrons as 
interacting via a meson exchange potential [3]. Hadronic, as opposed to quark, degrees of freedom 
should be appropriate for these situations. In this paper, we consider the S'-wave radial Schrodinger 
equation, 

h 2 d 2 



+ V(r) - E 



u(r)=0 , (1) 



2Mdr 2 

where M is the reduced mass of the two hadrons and u(0) = u(oo) = 0. 

In Section II, our focus will be on the form of the potential energy function V(r), especially 
as regards the contributions from multi-pion exchanges. In Section III, we carry out a numerical 
study of how varying the potential energy affects deuteron binding. Our conclusions and plans for 
future study are presented in Section IV. 



II. MULTI-PION EXCHANGE IN DEUTERON BINDING 

Although the deuteron state, and more generally the two-nucleon potential, are long-studied 
areas of research, they continue to attract theoretical attention. There has been much recent 
activity about how the nuclear potential behaves in the chiral limit of zero light quark mass [4-6]. 
This work has motivated our interest in deuteron binding. By design our description is simple, e.g. 
we omit tensor interactions and thus our 'deuteron' has no D-wave component [7]. We comment 
further on this aspect of our work in the Conclusion. 

Spatial potential energy functions for nonrelativistic calculations are obtained from particle 
exchange diagrams in quantum field theory. A nonrelativistic reduction of the exchange amplitude 
is followed by Fourier transformation from momentum to position space (see Eqs. (2)-(3) below). 
For example, photon exchange between an electron and a proton leads to a number of separate 
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effects (collectively the Breit-Fermi interaction [8]) starting with the long range Coulomb potential. 
In addition to generating long range potentials, such exchange amplitudes also generally give rise 
to short range or even local contributions. In this paper, we restrict our attention to only the long 
range part. 

A. Two-pion and Three-pion Exchange Potentials 

For convenience, we shall cast the potential energy in a dispersive form [9]. The strongest 
components of the two-nucleon potential arise from the isoscalar (I = 0) two-pion and three-pion 
exchange channels. A smaller isovector two-pion component must also be included. 

For isoscalar (1=0) potentials, the two-pion (scalar) exchange gives rise to an attractive inter- 
action while the three-pion (vector) interaction is repulsive. A momentum space version is given 
in terms of a corresponding spectral function pf~°\p), 



where i = S, V and /x| = (2m n ) 2 ,fi v = (3m w ) 2 . The spatial potential energie (r) corre- 

sponding to Eq. (2) is 



The above unsubtracted dispersive forms are taken from Refs. [6, 10], where behavior at relatively 
low energies (< 1 GeV) is studied. The issue of high energy behavior and subtractions is considered 
in Ref. [11]. However, the choice of a subtraction scale at sufficiently high energy is not expected 
to alter the basic conclusions of Refs. [6, 10], which are anchored by the constraints of chiral 
perturbation theory [12]. 

For the isoscalar vector channel, it is a good approximation to take p v ~°\p) oc 5(p, — m u ), 
where m u ~ 783 MeV. This is also appropriate for the isovector two-pion exchange channel, where 
now p v ~ l \p) (x 5(p — m p ), with m p ~ 770 MeV. In these two cases, the association between an 
exchanged vector meson of mass my and a spatial Yukawa potential Vy{r) oc e~ mvT is natural and 
obvious. The isoscalar two-pion exchange channel is more subtle. 

B. The Two-pion Isoscalar ('er(600)') Component of the Two-nucleon Potential 




(1=0) 



(2) 




(3) 



The nucleon-nucleon potential is known to have a strong attractive contribution of intermediate 
range. This arises from an exchange potential in the scalar, isoscalar channel, due to the exchange 



of a pair of pions in the I = S = channel. There is now a description of this which, although 
not explicitly given in terms of quark interactions, does nonetheless respect the underlying QCD 
dynamics. The approach is (i) to use chiral perturbation theory to describe the low energy part of 
the momentum-space potential, and (ii) to invoke two-pion rescattering via the Omnes function at 
higher energies in order to include the constraint of unitarity. 

The chiral component of p^~°^ turns out, as is generally the case in chiral perturbation theory, 
to depend on a small number of low energy constants whose numerical values are obtained phe- 
nomenologically. To the extent that these constants depend on any meson, it is actually p(770) 
that contributes. The chiral form of p^~°^ grows with energy and at some point requires modifica- 
tion in order to respect unitarity. Such a modification can be accomplished [6] by introducing the 
Omnes function, 



two ingredients of chiral perturbation theory and two-pion rescattering provide a sound theoretical 
basis for understanding the two-pion S = I = exchange potential. 

In the intermediate energy range 400 < E nw < 700 MeV, the S = I = two-pion phase shift 
does not pass through 90° and thus (unlike a; (783) and p(770)) has no obvious association with a 
resonant state. The dispersive spatial potential of Eq. (3) would seem to require a broad spectral 
function in this case. However, it can be shown (see especially Fig. 3(b) of Ref. [6]) that as a 
numerical recipe the effect of nonresonant S = I = two-pion scattering is reproduced by using 
V a (r) oc e~ m ° T with m a ~ 600 MeV. In our work, we adopt this approach and refer to it hereafter as 
the 'ct(600) exchange potential'. We stress that our use of '<j(600)' is simply a convenient shorthand 
and not meant to commit us to any specific model of this state [14]. 

III. BINDING THE DEUTERON 

Let us now consider the two-nucleon interaction in terms of four exchange potentials, 




(4) 



where 5^ °\s) is the S-wave two-pion phase shift in the 1 = channel. 1 In this manner, the 



y(NN) = y(NN) + y (NN) + yjNN) + y(NN) 



(5) 



1 It is also possible to proceed to a higher order of the chiral expansion [13]. 
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Each of these is taken to have the Yukawa form, 

„2 ra;r 

V (m (r) = Y,^ V ' (6) 

i 

where rji = +1 for repulsion and rfi = — 1 for attraction. All that remains is to fix the strength of 
each potential. 

For detailed discussions of potential energy fits to nuclear binding, see Refs. [10, 15]. We will 
assume the following values for the cr(600) and w(783) couplings [16], 

2 2 

and ~ 10.3 -> 12.9 . (7) 

47T 47T 

Our strategy will be to fit <7^ NN /47r to the experimental deuteron binding energy [17] 

BE dcut = 2.22457 MeV (with negligible error) (8) 

and see how it compares to the range in Eq. (7). We use this approach because a precise fit will 
allow us to study the effect on BEd eu t of modifying the potential energy function. 

The pion potential energy is obtained by starting from the field theory interaction with pseu- 
doscalar coupling 

£ttNN = ifnNNNj 5 T ■ 7tN (9) 

and reducing the pion exchange graph to its nonrelativistic limit. Upon keeping only the Yukawa 
dependence of Eq. (6), one obtains the isospin-dependent, spin-dependent potential energy 

*7rNN ( m K \ 

4vr \2M N J 

For the deuteron channel (1 = 0, S = 1) this reduces to 

K(r )(NN) = _£gN , (11) 



v^NN) (r) = ^/^y (2I 2_ 3) .H^_jij^ . fl()) 



with 
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^nn^nn(^) ^ 0.073 . (12) 

Observe that g^NN is tiny relative to g a NN and s^nn- 
For the p(770) interaction the lagrangian is 

£ p nn = 9 P mN^r ■ p^N . (13) 

where we take g 2 p /^ir = 2.0 [18]. The p(770) component of the potential energy is isospin dependent, 

(NN)(r) = V™ (2I 2 _ 3) eJW 
p 4-7T r 



1. The S=l, 1=0 (Deuteron) Channel 



Of the four component potentials, only the w(783) piece is repulsive in the 1 = channel. It is 
mainly the interplay of the large w(783) repulsion and cr(600) attraction that results in the deuteron 
binding. Although the coefficient of the u;(783) potential is larger than that of the attractive cr(600), 
the lighter mass of the <t(600) allows it to dominate over w(783) at intermediate distance scales. 
The w(783) provides the repulsive core. The potential energy which includes all four components 
is displayed in Fig. 1. When we use as input potential energies I/J NN \ Kj NN \ V^ NN ^ and fJ NN ^ to 
fit the value of <?^ NN / 47r to the deuteron binding energy, the value obtained for ^o-nn is 



a 2 

47T 



10.83 . 



(15) 



This is in accord with the range given above in Eq. (7). 
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FIG. 1: Deuteron Potential Energy. 



2. The S=0, 1=1 Channel 

The two-nucleon system has isospins / = 0, 1 and spins S = 0, 1, and if the two nucleons are in 
an S-wave, the I = S = and I = S = 1 configurations are ruled out by Fermi-Dirac statistics. 
The deuteron has / = 0, S = 1, which leaves finally the combination / = 1, S = 0. The I = 1, S = 
channel has a virtual bound state, as evidenced by the large scattering length a( 1 5'o) — 23.7 fm. 
Since <r(600) and w(783) are isoscalars, a potential energy with only these exchanges would imply 
equal binding in both the I = 0, S = 1 and / = 1, S = channels. However, the p(770) contribution 
is attractive in the deuteron channel but repulsive (at one third the strength) in the I = 1, S = 

a 



channel. Our numerical study obtains the desired result of no binding for the / = 1, S = channel 
upon using the <r(600) coupling strength of Eq. (15). 

A. Anthropic Implications 

Anthropics can be viewed as relating to the class of theories (including possibly string theory 
and chaotic inflation) in which spacetime has a domain structure, each domain with its own set 
of fundamental parameters (CKM mixing angles, lepton and quark masses, etc) [19]. Such a 
description is often referred to as the Multiverse [20, 21]. This theoretical possibility provides an 
interesting framework for viewing deep issues of physics, such as the cosmological constant [22], 
the Higgs vacuum expectation value [23] and the light quark masses [1]. 

Previous studies have revealed that even a modest variation of fundamental parameters can 
induce qualitative changes in physical systems [1, 23]. Works like these can help determine the 
'physical' windows in parameter space which enable Universes like our own to exist. This is in 
contrast with 'theoretical' parameter windows which can be associated with a given underlying 
theory, e.g. string dynamics. If the physical window is small compared to the theoretical window, 
then it is not unreasonable to assign a flat probability distribution across the physical window and 
thus to study the relative liklihood of that our Universe should exist within a given theoretical 
framework [20]. 

In this paper, we are concerned with the light quark masses. Refs. [1, 16] show that increasing 
iTiu + m-d decreases the aNN coupling whereas decreasing m u + rrid has the opposite effect. For 
example, it is estimated that taking — > produces a 40% increase in S^nn/^ 71 " while increasing 
m% to twice its physical value decreases S^nn/^ 71 " by 20%. The source of this sensitivity to m w 
(and ultimately to m u + lies mainly in the threshold contribution in the dispersive integrals of 
Eqs. (2), (3). By contrast, the pNN and wNN couplings are relatively unaffected by such changes 
in the light quark masses because their spectral functions are large only well above threshold. 

1. Unbinding of the Deuteron via Weakening of vJ NN ^: In our description, the two-nucleon 
potential is mainly the competition between the attractive V a and the repulsive V w . It is not 
hard to upset this balance and thus unbind the deuteron by varying some combination of the 
potential energies, e.g. by reducing V^ NN "* (keeping all other interactions fixed), increasing 
uJ NN -* (again keeping all other interactions fixed), etc. However, the former is the only 
variation associated with the mass values of the light quarks. To help keep track of how the 
deuteron binding energy responds to the variation in g a NN, we find it useful to define the 
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TABLE I: Dependence of deuteron binding energy on the criViV coupling as expressed directly (first two 
columns) or in terms of the ratios i?( cou P), i?( BE ) defined in Eq. (16). 

following ratios, 

[9.nn] = %^ [BE Acat ] = BEd ^ eV) . (16) 

From the numerical values in Table I, we see that the deuteron becomes unbound if V^ NN ' ) 
is reduced in strength by about 6%. This is to be compared with the 10% required to 
disassociate heavy nuclei [1]. 

2. Alternative Variations: There can, in principle, be other parameter variations. For the sake 
of completeness, we briefly consider a few of these. 

It is obvious that increasing will increase the deuteron binding energy. However, 

there is another consequence, which turns out to somewhat alter the physical world. We 
find that an increase of about 8% in V r J NN ' ) , with other potentials held fixed, suffices to first 
produce binding the / = 1,5 = NN channel. 

Varying the uj and p potentials would likewise have physical consequences, e.g. with 
Fj NN) fixed, an increase in strength of V^ NN ' ) by about 7% would result in unbinding the 
deuteron, This kind of variation is not associated in any obvious way with the fundamental 
parameters of the Standard Model and thus has no connections with our anthropic analysis. 

s 



IV. CONCLUSION 



Our aim in this paper has been to explore the dependence of deuteron binding on variations in 
m u + rrid in terms of a simple model involving just central potentials. The central potentials were 
constructed with care, by exploiting up-to-date results from studies of nuclear binding. This de- 
scription is admittedly in contrast with the traditional approach which uses the one-pion-exchange 
(OPE) tensor potential [7]. We feel that these two procedures are not incompatible. Consider a 
description containing the huge a and uj potentials along with the tensor OPE. If the a and uo 
effects cancel even more than in our calculation (which would involve a slight change in V a ), the 
deuteron binding would be given largely in terms of the tensor OPE. Even so, deuteron binding 
would still feel the changes induced by m u + rrid largely via the effect on V a (as described in this 
paper). We leave for future study a more ambitious description which includes all these elements 
in a Schrodinger equation context, plus insights from chiral perturbation theory. 

Thus, given the deuteron potential energies of Eq. (5), we have shown (cf Table I) that a 
slight reduction in V a results in unbinding the deuteron. Even a slight weakening of the attractive 
component results in a large effect. As a result, the deuteron is a sensitive system for studying 
how physics would respond to changes in m u + m^. The crucial links in the chain of logic which 
connects these are (i) the chiral relation, = B§{rn u + m^), which relates the pion mass to the 
light quark masses and (ii) the dispersive formula Eq. (2) which relates V a to the pion mass [10]. 
The 6% decrease in V a needed to unbind deuterons is less than the estimated 10% needed to unbind 
heavy nuclei [1]. 

It is abundantly clear that if conditions forbade atoms or nuclei from existing, then there would 
be no Universe as we know it. The absence of the deuteron, although itself just a solitary quantum 
particle, would likewise be dramatic. The deuteron first appears in Cosmology as part of Big Bang 
Nucleosynthesis (BBN) in the early Universe, 



Before deuteron formation, free neutrons are copious (about one neutron per five protons), but after 
formation free neutrons are absent (the number of free neutrons equals the number of deuterons 
at about 200 seconds after the Big Bang). Once deuterons are formed, they take part in a number 
of additional BBN reactions, such as fusion with protons and neutrons, 



p+ n # D + 7 



(17) 



D+p ^ 3 He + 7 D + n ^ 3 H + 7 , 



(18) 



Q 



or with each other, 

D + D # 3 H+p D + D # 3 He + n . (19) 
These lead to the production of 4 He, 

3 H + p # 4 He + 7 3 He + n # 4 He + 7 

(20) 

3 H + L> # 4 He + n 3 He + L> # 4 He + p 

Although trace amounts of heavier nuclei are formed, the above relations display the basic yield of 
BBN once the deuteron is formed. It is clear that, although just a 'solitary' particle, the deuteron 
occupies a central role in BBN. 2 

There are some additional issues, not yet mentioned, which we leave for further study. As 
regards BBN, it would be interesting to study more carefully what happens under conditions 
where the deuteron becomes slightly unbound. In a different direction, by using a modification of 
the meson exchange approach described in this paper, we can consider the possibility of molecular 
structures in the meson sector. Of particular interest are the states do (980) and /o(980), whose 
nearness to the KK threshold has suggested a molecular interpretation [24]. We will report on 
this study elsewhere [25]. 
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